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I. INTRODUCTION 



Lattice QCD is a unique tool to compute the mass spectrum, leptonic decay constants 
and hadronic matrix elements of local operators non-perturbatively from first principles. 
Renormalization of lattice operators is an essential ingredient needed to deduce physical 
results from numerical simulations. In this paper we study the renormalization properties 
of composite bilinear operators with the overlap quark action. 

In principle, renormalization of quark bilinears can be computed by lattice perturbation 
theory. However, it is generally difficult to go beyond one loop in such calculations. To 
overcome these difficulties, Martinelli et al. flj have proposed a promising non-perturbative 
renormalization procedure. The procedure allows a full non-perturbative computation of 
;he matrix elements of composite operators in the Regularization Independent (RI) scheme 

m 

. The matching between the RI scheme and MS, which is intrinsically perturbative, is 
computed using only the well behaved continuum perturbation theory. 

This method has been shown to be quite successful in reproducing results obtained by 
other methods, such as chiral Ward Identities P]. The method has also been successfuUv 
applied to determine renormalization coefficients for various operators using the Wilson |4 , 
^ fl El, staggered j^, domain- wall 0, chirally improved 10 1, and overlap fermions [h, 
12| |. The purpose of the current work is to study the application of this non-perturbative 
renormalization procedure to the renormalization of the quark field and the flavor non-singlet 
fermion bilinear operators, and also to study their quark mass dependence for the case of 
overlap fermions. 



Neuberger's overlap fermion 



symmetry on the lattice 13|, 



IM, with flnite cut-off. As a consequence, many chiral- 
symmetry relations and the quark propagator 17] preserve the same structure as in 



14| is shown to have correct anomaly and exact chiral 



the continuum. The use of the overlap action entails many theoretical advantages |18|: it 
has no additive mass renormalization, there are no 0{a) artifacts, and it has good scaling 
behavior with small O(a^) and O(m^a^) errors In addition, it forbids mixing among 

operators of different chirality and, therefore, can be very helpful in computing weak matrix 
elements. 

The outline of this paper is as follows. In Sec. m we review the non-perturbative method 
(NPM) proposed in Ref. jl| and introduce the notation used in the remainder of this work. In 
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Sec. mil we briefly describe the overlap fermion formalism. We present the numerical results 
for the renormalization constants as well as the Renormalization Group (RG) analysis of the 
quark bilinear in Sec. IIVI and Sec. In Sec. IVIl we examine the flnite m behavior of the 
renormalization factors of the quark bilinear operators. We complete our discussion with 
our conclusions in Sec. IVlll 



II. NON-PERTURBATIVE RENORMALIZATION METHOD 

In this section, we review the nonperturbative renormalization method of Ref . Q , which 
we will use to compute the renormalization constants of quark bilinears in this paper. The 
method imposes renormalization conditions non-perturbatively, directly on quark and gluon 
Green's functions in Landau gauge. 

We start by considering the deflnition for the momentum space quark propagator. Let 
S (x, 0) be the quark propagator on the gauge-flxed conflguration from a source to all 
space-time points x. The momentum space propagator is deflned as the discrete Fourier 
transform over the sink positions 

S (p, 0) = ^ exp (-zp'^" ■ x) S {x, 0) , (1) 

X 

where p'*^** is the dimensionless lattice momenta. 

In our case, we use the periodic boundary condition in spatial directions and the anti- 
periodic boundary condition in time direction. We have then the dimensionful momenta 

271 2n 
p. = -— (n, - iV,/2), and pt = —{nt - 1/2 - Nt/2) , (2) 

for an x Nt lattice. 

We also deflne the square of the absolute momentum as the Euclidean inner product of 
the momenta defined in Eq. ((21) 

where we use convention that p is dimensionful and p)^" is dimensionless. 
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A. Three Point Function 

Consider the flavor non-singlet fermion bilinear operator 

Or(x) = i,{x)Ti,, (4) 

where F is the Dirac gamma matrix 

r G {1,7m,75,7m75,c^mv} , (5) 

and the corresponding notation will be {S, V, P, A, T} respectively. The flavor index is 
suppressed. The connected three point function with an operator insertion at position 
between the quark flelds at x and y is given by 

Go(x,0,i/) = (^(a;)Or(0)^(?/)) = 0)15(0, y)), (6) 

where 5'(0, y) is the quark propagator from y to 0. It is the inverse of the Dirac opera- 

averaging over all 

gauge conflgurations, i.e., 

(5(x,0)), (7) 

is it translational invariant. 

Using 75 hermiticity, the Fourier transform of the three-point function is given by 



= {(^j d^xSix,0)e-^P-^^T(^J dV(0,i/)e*P'-^^), 
= {S{p, 0)F75 (1 d'yS^iy, 0)e^^'-^^ 75), (8) 
where the f refers only to the color and spin indices. This can be written as 

Go{pa,p'a) = (5(p,0)F {i,S\p' ,0^,)), (9) 
From this, one can deflne the vertex function as the amputated three-point function 

Ao{pa,p'a) = S{pa)^^Go{pa',p'a')S{p'a)^^ , (10) 

where 

S{pa) = {S{p,0)), (11) 



which is translational invariant and is a 12 x 12 matrix in color-spin space. 
Finally, a projected vertex function is defined 



ro(pa) = —l-Tr (Ao{pa,pa)Po) , (12) 
Tr(P^) V / 

where Pq = F is the corresponding projection operator. 

B. RI-MOM Renormalization Condition 

The renormalized operator 0(/x) is related to the bare operator 

0(/i) = Zo{f^a,g{a))0{a) , (13) 

and the renormalization condition is imposed on the three-point vertex function Fo(pa) at 
a scale = fi^ as 

r ^ M Zo{fia,g{a)) 

Fo,ren(pa) p2=^2 = — -—T o{pa)\p2=^2 = 1 (14) 

Z^{fj,a,g{a)) 

n 

to make it agree with the tree-level value of unity |1|] . Here is the field or wave- function 
renormalization 

l/jren = ^^V- (15) 

In order to alleviate the non-perturbative effects from the spontaneous chiral symmetry 
breaking, high virtuality with fi ^ ^qcd is required. On the other hand, to avoid the 
discretization errors, one would need /i <^ 1/a. So, for the RI-MOM procedure to be a valid 
and practical renormalization scheme, there should exists a window in the renormalization 
scale /i, i.e. Aqcd ^ ^ 1/a. In this work, the renormalization constants are extracted 
from {pa^ > 1 which, in principle, should have large discretization errors. But, as we will 
see later, we do not see large discretization errors. This is also the case observed in previous 
studies |l|,y,y,bj. 

In practice, one usually matches the results to the perturbative scheme, e.g. MS scheme, 
in order to compare with experimental quantities. We will discuss perturbative matching 
in Sec. IIVI In general, the vertex function Toipa) may have intrinsically non-perturbative 
contributions, e.g. from the Goldstone boson propagator, which are not included in pertur- 
bative calculations. To this end, we either go to large enough momentum such that the 
non-perturbative effects are suppressed or somehow remove them from the data. 



There are several ways to obtain the renormahzation constant Zq for the operator O. 
For the Z^ in Eq. (fT^ . one could use a known ratio involving Z^ to have it eliminated in 
Eq. (|14p. The first way is to extract Z^ from the vertex function of the conserved vector or 
axial- vector current. For example, if one uses the conserved vector current, then Zy^-, = 1. 
From the renormahzation condition 



ryc(P«)lp2=M2 



one then obtains 



(16) 



(17) 



One can also extract Z^ directly from the quark propagator. From Ward Identity (WI), 
it follows [l| 



dS{pa)-^ 



(18) 



To avoid derivatives with respect to a discrete variable, it is suggested [1| to use 



-«Tr^f;^i7^(p/,a)5(pa) ^ 
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(19) 



which, in the Landau gauge, differs from Z^ by a finite term of order aj. [20]. The matching 
coefficient can be computed using continuum perturbation theory, and up to order 

^ = 1 - y 
Z^ 



(47r 



2 AA^ -h . . . 



(20) 



In the Landau gauge 



a(2) = (3 + 22iV2 - 4N,nf) 



(21) 



^ 16 

where is the number of colors and Uf the number of dynamical quarks. However, as 
pointed out in Ref. [sl, due to the ambiguity as to how the discrete lattice momentum p is 
defined, this method will introduce roughly 10% - 20% uncertainty in determining Z^. 
The third way is to first calculate the renormahzation constant Za from the axial Ward 



identity 
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m 



axial current 



with a local current, and then use the renormahzation condition for the 
Za- 



■'ill 



-TAipa)l 



1, 



(22) 
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to eliminate the unknown from Eq. ()14|) . Combining Eqs. ()14|) . ()22p and the from 
the Ward identity, one obtains 

Zo = Z^Jf^. (23) 

This way, other renormahzation constants, such as Z5, Zp, Zy, and can be obtained and 
the identity relations Zs = Zp and Zy = Za due to chiral symmetry can be checked. Note 
here, the local axial-vector current is finite. Thus, Za is independent of scale, but depends 
on the lattice spacing a. 

In this work, we shall adopt the third approach as mentioned above. It is known that Za 
as determined from the Ward identity has a small statistical error at the level of 0.2% 21] and 
will not contribute much to the overall error. When we study the quark mass dependence, 
we shall use the wave-function renormahzation from the quark propagator to obtain its 
dependence for the renormahzation factors of the quark bilinear operators. 

III. OVERLAP FERMION 

The massless overlap-Dirac operator in lattice units is Q| 

D(0) = [l+75e(i/)] , (24) 

where €{H) is the matrix sign function of an Hermitian operator H. e{H) depends on the 
background gauge field and has eigenvalues ±1. Any such D is easily seen to satisfy the 
Ginsparg-Wilson relation [2^ 

{75,D} = D75^. (25) 

For the topological sector with no zero modes, it follows easily that {■y^, D~^{0)} = 75 and 
by defining Z)^^(O) = [D^^{0) — 1/2] we see that it anticommutes with 75 

{75,£>"'(0)} = 0. (26) 

The standard choice of e{H){x,y) is e{H) = H]y/\Hw\ = / {hI^^HwY^'^ , where 
Hw{x,y) = -y5Dw{x,y) is the Hermitian Wilson-Dirac operator. Dw is the usual Wilson- 
Dirac operator on the lattice. However, in the overlap formalism the Wilson mass parameter 
p needs to be negative in order to generate zero modes. 



In the present work, we use the standard Wilson-Dirac operator, which can be written 

as 



Dw{x,y) 



(27) 



The negative Wilson mass —p is then related to k by 

1 



K = 



(28) 



2(-p) + 8 ■ 

p is chosen such that k > Kc and p < 2r, and hence there is no species doubling and there 
can be zero modes for the massless quark. We take r = 1 and p = 1.368 (which corresponds 
to K, = 0.19) in our numerical simulation. 

It is shown that the flavor non-singlet scalar, pseudo-scalar j^], vector, and axial j^jQ] 
bilinears in the form tjjKT(l — ^D)ip {K is the kernel which includes 7 matrices and T 
is the flavor SU{Nf) matrix) transform covariantly under the global chiral transformation 
Sip = T75(l — D/2)iIj as in the continuum. The 1 — factor is also understood as the 
lattice regulator which projects out the unphysical real eigenmodes of at A = 2. For the 
massive case, the fermion action is deflned as ippDif) + matp{l — \D)il) so that the tree-level 
wave-function renormalization of the quark propagator is unity. In this case, the Dirac 
operator can be written as 



D{m) = pD + ma{l - -D) 



(29) 



With the ip fleld in the operators and the interpolation flelds for hadrons replaced by the 
lattice regulated fleld ^/^ = (1 — \D)ip, the regulator factor will be associated with the quark 
propagator in the combination (1 — |i5)D(m)~^ in Green's functions, leading to an effective 
quark propagator jl?! 

1 



1 



-D)D(m) 



m] 



(30) 



where the operator Dc = pD/{l — ^D) is chirally symmetric in the continuum sense, i.e. 
{75, Dc} = j2^,l3|; but, unlike D, it is non-local. Thus, the effective quark propagator 



in Eq. (jHlHl turns out to have the same form as in t 



le continuum, i.e. a chirally symmetric 



Dc plus a mass term in the inverse propagator 



relation of the pseudoscalar and vector mesons, it is learned 



. By studying the dispersion 
1711 that the (ma)'^ errors are 



much smaller than those of the Wilson fermion, making it a viable option for studying 



the heavy-light systems. Furthermore, it affords a non-perturbative renormahzation of the 
heavy-hght decay constant via the chiral Ward identity and the unequal mass Gell-Mann- 
Oakes-Renner relation jl?!. The preliminary study of the charmonium spectrum with the 
quenched overlap fermion seems encouraging as far as the hyperfine splitting and the S-wave 
to P-wave charmonium splittings are concerned jsol - 



IV. NUMERICAL RESULTS 

In this paper we work on a 16^ x 28 lattice with lattice spacing, a=0.20fm, as determined 
from the pion 

gauge action through the pseudo-heat-bath algorithm. A total of 80 configurations are used. 
The lattice parameters are summarized in Table HI 



TABLE I: Lattice parameters. 



Action 


Vokime 


ATherm 




(3 


a (fm) 


Physical Volume (fm'*) 


Iwasaki 


16^ X 28 


10000 


5000 


2.264 


0.200 


3.2^ X 5.60 



Recently it has been speculated that the overlap operator with coarse lattice spacing 
of 0.2 fm, such as in Ref. 0] and this study, might have a range as large as 4 lattice units. 
Thus, the calculations rnight be afflicted by unphysical degrees of freedom as light as 0.25 
GeV. It has been shown j33| by direct calculations at lattice spacings of 0.2 fm, 0.17 fm, and 
0.13 fm that this is not the case. The range of the overlap operator as defined in Ref. 3^ 
is about 1 lattice unit (in Euclidean distance or 2 units of 'taxi driver' distance) for each of 
the above lattice spacings. Therefore, the range of the overlap Dirac operator scales to zero 
in physical units in the continuum limit and our 16^ x 28 lattice, with spacing of 0.2 fm, 
used in this study is in the scaling range. 

The gauge field configurations are gauge fixed to the Landau gauge using a Conjugate 
Gradient Fourier Acceleration js^ algorithm with an accuracy of = ^ |9^^^(x)|^ < 10^^^. 
We use an improved gauge-fixing scheme to minimize gauge-fixing discretization errors. 
Since gauge-fixing is involved for the external quark state, there are concerns about the 
effects of Gribov copies on the numerical results of the renormahzation procedure. A study 
of two different realization of the Landau gauge and a covariant gauge shows that the 
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renormalization constants from these gauge-fixings differ by less than the statistical errors 
of about 1 - 1.5% level jjsTj. We shall thus assume that that the potential uncertainty due 
to Gribov copies is at this 1 - 1.5% level which is comparable to our statistical errors. 

Our numerical calculation begins with an evaluation of the inverse of D{m) which is 
defined in Eq. (|29|) . We use a 14th-order Zolotarev approximation [38] to the matrix sign 
function e{Hw)- In the selected window of x G [0.031,2.5] of e(x), the approximation is 



better than 3.3 x 10 



-10 



We then calculate Eq. (jHUj) for each configuration by using 
multi-mass Conjugate Gradient method for both the inner and outer loops. The detailed 
numerical description is given in Ref. 0- 

In the calculations, k = 0.19 was used, which 
corresponds to p = 1.368. We calculate 15 quark masses by using a shifted version of the 
Conjugate Gradient solver jiol . Iigj ). The bare quark masses ma are chosen to be ma = 
0.02100, 0.03033, 0.04433, 0.06417, 0.07583, 0.08983, 0.10850, 0.12950, 0.15633, 0.18783, 
0.22633, 0.26833, 0.32200, 0.40000, and 0.60000. With the scale determined by they 
correspond to pion masses 212(7), 247(6), 290(6), 342(6), 370(7), 400(7), 438(7), 478(8), 
524(8), 575(9), 633(10), 692(11), 764(12), 862(13), 1092(17) MeV respectively 
In the following, we give the steps for the numerical calculation: 

• After we calculate the quark propagators in coordinate space for each configuration, 
we use the Landau gauge transformation matrix to rotate the quark propagators to 
the Landau gauge. The discrete Fourier transform is then used to calculate the quark 
propagators in momentum space. 

• Next, we calculate the five projected vertex functions To{pa) defined in Eq. (fT^ . 
where we have used the effective quark propagator in Eq. (jHUj) for the calculation. 
By definition, they are the ratios of renormalization constants at the chiral limit (i.e. 
Z^{fia, g{a))/ ZoifJ'a, g{a)) from Eq. (IHj)) in the RI scheme at scale /i^ = p"^. They are 
in general dependent on (pa)^ which comes from two sources. One is from the running 
of the renormalization constants in the RI scheme; the other is from the possible (pa)^ 
error. 

• We decouple the two above mentioned scale dependencies of the calculated Toipa) 
by first dividing out its perturbative running in the RI scheme. Ideally, this should 
take care of the scale dependence, since we have taken the scale to infinity. However, 
due to the (pa)^ error on the lattice, there can still be some (pa)^ dependence in the 
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To{pa) after undoing the perturbative running. Following Ref. [9f, we shall attribute 
the remaining scale dependence to the (pa)^ error and will use the simple linear fit 
to remove it. This will be discussed in Section For the scalar and pseudo-scalar 
vertex functions, there is an additional complication due to the presence of quark mass 
poles 0]. We shall remove them first and then extrapolate to the chiral limit. Finally, 
we can check the expected relations Za = Zy and Zs = Zp . 

• In order to compare results with experiments, one frequently quotes the scale- 
dependent results in the MS scheme at certain scale. So the final step is to per- 
turbatively match the results from the RI scheme to the MS scheme at /i = 2 GeV for 
Zgj Zp Zj'j and Z^. 



A. Axial and vector currents 



Let us consider first the vector and axial vector currents. Since each obeys a Ward 
Identity .39], their renormalization constants are finite. In Fig. ^ we show the vertex 
functions and Ty from Eq. (fT^ for different bare quark masses as a function of the 
lattice momentum {pa)"^. We find that they are weakly dependent on the mass, and almost 
scale independent after (pa)^ > 2.0, which corresponds to p > 1.4 GeV. 

In the RI scheme, Eq. fll4j) implies that in the chiral limit 

]im r^/y(pa)|p2=^2 = Z^{ij,a,g{a))/ZA/v, (31) 

and one expects that Za = Zy for the overlap fermion, but this is true only for large 
momenta p. At low momenta, F^ and Fy may differ due to the effects of spontaneous chiral 
symmetry breaking Q]. 

Following Ref. we show in Fig. |2lthe quantities F^ — Fy and |(FA + Fy), after linearly 
extrapolating with respect to ma to the chiral limit. We can observe from the upper panel 
of the figure that there is no effect of spontaneous chiral symmetry breaking at moderate 
and high momenta, where F^ — Fy tends to zero. In fact, the percentage error between 
F^ and Fy is less than 1% at (ap)^ = 2.5 (top panel) and is comparable to the statistical 
error. This deviation decreases further for higher momenta. On the other hand, the effects 
of spontaneous chiral symmetry breaking are clearly visible at low momenta where F^ and 
Fy differ. In the lower panel of Fig. El we plot |(F^ + Fy) against {pay. Since Zy and Za 
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FIG. 1: The projected vertex function T defined in Eq. 1)12(1 for the vector and 
axial- vector currents with different bare quark masses. 

are scale independent, the slight (pa)^ dependence that one observes in |(r^ + ry) in Fig. |21 
for (pa)^ > 2.5 reflects the scale dependence of and the lattice (pa)'^ error. 



B. Pseudo-scalar and scalar densities 

The pseudo-scalar and scalar densities differ from the axial and vector currents in the 
sense that their renormalization constants are not scale independent. The scalar and pseudo- 
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FIG. 2: The upper panel is {Ta — Tv) versus (pa)^. We see that Za = Zy 
is valid above moderate (pa)^ ~ 2.5. Here we linear extrapolate to the chiral 
limit (m = 0) with respect to ma. The lower panel is ^ (Ta + Ty) versus 
(pa) 2. 



scalar densities with the form ?/'(l — D/2)ip and ijj'y^il — D/2)iIj transform under the lattice 



chiral transformation as in the continuum Il5 



4G 



4l|. From the Ward identities, one obtains 
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the relations 



Zs = Zp (32) 
Zm = ^. (33) 

Thus, the quantities Zs/Zp, ZgZm and Z^Zp are expected to be scale independent. 

For the case of the pseudo-scalar and scalar renormalization, there is a complication due 
to the presence of the quark condensate in the inverse quark propagator. Using the axial 
Ward identity from the quark propagator, one has Q 

mVp{j>,V) = YTr{S-\v))- (34) 



It is known that due to the spontaneous chiral symmetry breaking, the 



;race of the inverse 



quark propagator picks up a contribution from the quark condensate j4^ . At large j^l , 
it is given by 



-Tt{S-\v)) = m- {qq)^ + 0(1//). (35) 



^Tr{S \p)) = m- (gg)^' 

n 

from first order perturbation j42]. This implies that the renormalized Tr{Sj.J^{p)) should be 



^TriS;M) = mren - + 0{l/p') . (36) 

12 



In the study of lattice artifacts of the Wilson fermion 



43|, it is shown that there are three 



terms which mix at 0{a) to give an improved and renormalized quark field, 

qren = ^^^(1 + bqma){l + ac'gip + rriren) + acNGi ^}go, (37) 

where P may appear due to gauge fixing. It is found in the study of the order 0(a) 
error of Wilson fermion, that is large. Combining Eqs. (j36j) and (j37p . one has Q] 

J^Tr■{Sf^l^{pa)) = ... + ZmZ^{ma+mresa)-CiZ^ ^^^^^ +2{cngi- Cq){paf + 0{l /p'^), (38) 

where terms of OimcNci) are neglected. Thus, in this case, j^Tr{Si^l^{pa)) diverges for 
large pa. 

On the other hand, it is learned from the domain wall fermion study on a 16'^ x 32 x 16 
lattice with Wilson gauge action ai j3 = 6.0, the explicit chiral symmetry breaking effect 
from the CNGi~c'q term is negligible for moderately large values of (pa)"^ and that the residual 
mass is small. Since the explicit symmetry breaking is controlled to a level < 10~^ with 
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FIG. 3: A plot of nrr(S';;^\(pa)) versus (pa)^ for different bare quark mass ma on 16^ X 28 lattice, 
showing that j2Tr{S^J^^{pa)) approaches a constant value at large (pa)^. 



the Zoloterav approximation of the sign function [31| for the overlap fermion, we expect the 
cngi ~ c'q term to be negligibly small. It was already shown that the residual quark mass 
due to the numerical approximation of the overlap operator is neglig ible Q. Here, we shall 
verify the expectation that the Cngi — c'q term is indeed negligible. 

As shown in Fig.lHl j2Tr(Si^l^{pa)) for several quark masses tend to constant values after 
{pay > 4 . This indicates that there is no discernible contamination due to the explicit 
chiral symmetry term 2{cngi — d^{paf' which grows as {paY- Fig. E] shows the same at the 
chiral limit which is obtained from linear extrapolation in ma. In this case, j^Tr{S^l^{pa)) 
tends to zero at large {pa)"^ as expected from Eq. ()38|1 with no residual mass. Plotted in 
Fig. El is the slope of j^Tr{Si^l^{pa)) with respect to the quark mass ma. It is expected to 
be ZjnZ^ at large {pa^ from Eq. ()38|1 . These results are very similar to those of the domain 
wall fermion Q]. 

From Eqs. and (|HHj). one finds that 

Zj 



Tp{pa,pa) 



c,z,^-P- + o{W) 

ma[pa)'^ 



(39) 



Since the quark condensate {qq) has a contribution of 4^ from the zero modes due to 



mV 



the topological charge Q, it is expected that rp{pa,pa) has 1/m^ and 1/m singularities as 
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FIG. 4: The value of j2'^''' {Si J^ipa)) extrapolated to m = vs (pa)^ by a simple linear extrapo- 
lation. At large (pa)"^, the extrapolated value is zero within errors. 

m ^ 0. Thus, it is suggested jo^ to fit rp{pa,pa) with the functional form involving pole 
terms. As illustrated in Fig. IHl the singular behavior in m is quite visible. It is suggested in 
Ref. Pi to fit the Fp with a double and single pole form for each pa 



Tp{pa,pa) 



Cl,P C2,P , / \2 

+ r + C3,p + CApiam) 



(40) 



[am)^ (am) 

This is appropriate for the case studied in Ref. where the lattice volume is relatively small 
(the space-time volume is ~ 10 fm'*) so that the zero mode contribution is substantial and 
the quark mass is relatively heavy so that the quenched chiral log is not significant. In our 
case, the space-time volume at 184 fm^ is much larger. As such, the zero mode contribution 
is expected to be smaller. In our study of the quenched chiral log in the pion mass with 
the same lattice, it is found that the pion mass is basically the same when calculated from 
either the (PP), {A4P), {A4A4), or {PP — SS) correlators, indicating that the zero mode 
effects are small and negligible within statistical errors, even for the smallest pion mass at 
~ 180 MeV. Our lowest pion mass is 212 MeV in the current study. On the other hand, 
the quenched chiral log is quite prominent with pion mass less than 400 MeV. Therefore, we 
believe that the more appropriate approach is to relate the quark condensate in Eq. (j39j) to 
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FIG. 5: A plot of the slope of j2'^''~ i^iattiP'^)) with respect to the quark mass as a function of 
(pa)^. It is expected to be ZmZ^ at large (pa)^ from Eq. (|38|) . 



the pion mass through the Gell-Mann-Oakes-Renner relation 



2m 



(41) 



and use the power form for the pion mass where the leading log is re-summed through the 
cactus diagrams [4^, i.e. 



m^(ma) = A (ma) ^+'' + B{maf . 



(42) 



Here 5 is the quenched chiral log parameter. 

Neglecting higher order terms and expanding ^ /7r(0) + c(ma), we can approximate 
Eq. dnni) as 



rp|^^^(ap,ma) ~ AA{ma)^+i> ^ + ap{ma)^+s ^ > + A2 + Cp{ma ) + Dp{may , (43) 



where 



Ai 
ap 

A2 



2 (pa) 2 
c 



CiZ 



MO) ' 

2 (pa) 2 



(44) 



17 




FIG. 6: Fp versus (pa)^ with different masses. Here, one can clearly see the 
strong divergent behavior for small quark masses. 
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FIG. 7: Ts versus (pa)^ with different masses. The chiral log behavior is also 
quite visible. 

We plot rp{pa,pa) in Fig. Elas a function of (pa)^ for several ma. It is clear that at small 
(pa)^, it has a singular behavior for small ma which is presumably due to the divergent 
terms associated with the quenched chiral 6. 
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From the vector Ward identity, one has the relation 

1 dTr[S{p)-'] 



12 dm 

From Eq. (jHH|l . the above can be approximated with 



(45) 



r5(pa,pa) = - + ^^-— + ... (46) 

for large [pa)"^. 

Similar to the pseudo-scalar case, substituting the pion mass in Eq. fl42|l and the Gell- 
Mann-Oakes-Renner relation in Eq. (j41|) to Eq. (j45|) and neglecting the higher order terms 
beyond (ma)^, we obtain 

r5,latt(«P'"^«) = Y^M^"' + «SY^("^a)^"'| 

+A2 + Cs{ma'^) + Ds{maf, (47) 

At first glance, it appears that the quantity Z^Z^ in A2 is not separable from the other 
term in Eq. ()44|1 . However, we should note that we know the values of A,B and /7r(0) in 



Eq. (j44j) from an earlier study of the pion mass and decay constant [3l|. After fitting Ai 
and A2 in Eqs. (gSI) and (jUj), one can compare the first term in A2 with Ai to obtain ZmZ^. 
As we shall see later, it turns out the first term in A2 is 0(10^^) times smaller than Z^Z^. 

In order to obtain Z^Z^ and assess its finite m behavior, we first subtract the divergent 
terms in Eqs. (jlHj) and (jTfj) and then fit the subtracted vertex functions linear and quadratic 
in m, i.e. with ma^ and m^a^ terms. In the following sub-section, we shall detail our fitting 



methods and give the results. We should note in passing that it 



non-degenerate quark masses to avoid the divergence in Fp 



11 



las been suggested to use 
However, due to the 



complication of the quenched chiral log, it is not applicable here. 



1. Fitting 

We adopt the fitting procedure used to fit the chiral logs in the pion mass j,31l] and the 
Roper resonance in the nucleon correlator with priors. From the chiral log fit of the 
pion mass j^^, we obtain 6 to be in the range of 0.20 - 0.15 when the maximum pion mass 
for the fitting range is set to be ~ 500 MeV - 900 MeV. Since we are fitting Fp and F5 in 
the similar quark mass range, we put a weak constraint on the value of 6 with 6 = 0.18(5) 
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TABLE II: Fitted parameters corresponding to Eqs. (|43() and H47() 



Channel 


Ai ap{as) S 


C 


D 


xVdof 


P 


0.008(1) 0.08(9) 0.163(18) 1.305(15) 


-0.024(6) 


-0.022(5) 


1.29 


S 


0.011(5) 0.14(42) 0.171(24) 1.302(17) 


-0.079(8) 


-0.020(5) 


1.47 



which covers the range of 5 in the fit of the quenched chiral log in the pion mass. Data 
corresponding to the few lowest masses are first fitted with A\p{A\s)^ ap{as) and A2p{A2s) 
and then these parameters are constrained with those fitted values to fit the whole range of 
the masses with the forms in Eqs. (|43p and (j47|) . It is observed that these vertex functions 
are highly correlated between different masses and that the correlation increases with higher 
momentum. This, we believe, is due to the fact that the quark masses that we are interested 
in are all much smaller than the external momentum of pa = 4.145 which corresponds to 
/i = 2 GeV that we will use to eventually match to the MS scheme at this scale. In this 
sense, the high correlation is a generic feature that this non-perturbative renormalization 
procedure faces for hght quarks. 



We show in Fig. IHlthe vertex functions rp{pa,ma) and rs{pa,ma) at pa = 4.145 as a 
function of ma, together with the fitted curves based on Eqs. PHjl and (jUj). The fitted 
parameters are given in Table ITTl We see from the fitted parameters Ai,ap{as),S and A2 
from Fp agree with those fitted from F^ respectively, as we expected. This supports our 
supposition that the singular behaviors in both the Fp(pa,ma) and rs{pa,ma) are due to 
the quenched chiral log in (qq). We also tried to fit the pseudoscalar vertex function in the 
form in Eq. pUj) and found that it does not fit well — the is too large. 

We note that the fitted central values of 6 tend to be on the low side compared to our 
previous fit of the pion mass which gives 6 = 0.20(3) jsil and is in agreement with 6 ~ 0.23 
as deduced from the topological susceptibility calculation with the overlap operator j^|. 
We think the reason is that the smallest quark mass in the present fit which corresponds to 



m-TT ~ 250 MeV is larger than that in the previous fitting o: 
to ~ 180 MeV. According to the detailed study 



the pion mass which corresponds 
of the quenched chiral log as a 



function of the fitting range of quark masses, this behavior of a smaller 6 for a higher quark 
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FIG. 8: Vertex functions for scalar (top) and pseudo-scalar (bottom) channels at the momentum 
corresponding to /u = 2 GeV. 

mass range is to be expected. 

When we take the value of A = 1.3 and B = 1.1 from our pion mass chiral log fit (Fig. 12 



in Ref. 



31[) in the relevant mass range, it follows from Eq. ()44|) that 



2(paf A 

which is two orders of magnitude smaller than A2. We shall subtract this contribution from 
A2 in Eq. (j44j) to obtain Z^Z^ which changes its value by about half a a which is not 
significant. 

To eventually obtain Zs and Zp and their respective finite m dependence, we define the 
subtracted Fp and F5 by taking out the divergent terms in Eqs. and and the first 
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term in A2 in Eq. ()44|) on each Jackknife sample (J) as 
r5,"*'-^(ma) = r^(ma) - Af;|(ma)^ 



+ aj (ma) 



-1 



A 



[ma 



l+S, 



(49) 
-(50) 
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FIG. 9: Z^/Zp and Z^jZg and the ratio Zp/Zs as a function of (pa)^. 

From Eqs. PHj) . (jUj), (jlTj) . and (jl^ . we see that the ratios of Z^/Zp and Z^jZg are the 
subtracted vertices Fp and F5 at the massless hmit for each -pa in the RI scheme, i.e. 



F^"^(pa,ma = 0). 



(51) 



^p,5(pa) 

We plot in Fig. El Z^jZp, Z^/Zs, and the ratio ZsjZp as a function of (pa)^. We see 
that for (pa)^ > 3, the ratio goes to unity which is a confirmation that our fitting procedure 
does not spoil the expected chiral property Zp = Zs for the overlap fermion. 



C. The tensor current 



In Fig. EH we show Ft versus (pa)"^ with different masses. We can see that at moderate 
to large (pa)^, F-r is not sensitive to the quark masses. The chiral limit value is obtained 
by a linear plus quadratic fitting as for the case of vector and axial vector currents and the 
results will be presented in the next section. 
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V. RUNNING OF THE RENORMALIZATION CONSTANTS 



In general, one can choose to define the renormahzation conditions for different p and p' 
in Eq. (jiup . But the virtuahty of the quark states must be much larger than Aqcd- This is 
so because in order to obtain physical results at certain scale (e.g. momentum or mass), one 
needs to combine the matrix element of the renormalized operator 0{^) with the Wilson 
coefficient function. The latter is usually computed in continuum perturbation theory by 
expanding in at a scale of order of fi. Thus, for the validity of perturbation calculation, 
/i must be large. On the other hand, one would like to have /i <^ 1/a in order to have 
smaller 0{a) effects. When spontaneous symmetry breaking takes place, as is in QCD, a 
large /i may not be enough due to the presence of the Goldstone boson. For example, at 
low momentum transfer q = p' — p, the Green's function can have a Goldstone boson pole 
like However for fermions (but not for a scalar particle), this contribution will be 

= smaller than the perturbative contribution even when = Q]. Thus, it 
is desirable to have a window Aqcd <^ yU ^ 1/a so that both the non-perturbative effects 
and the lattice artifacts are small. In practice, one finds that the renormahzation procedure 
prescribed here works well for (/xa)^ as large as 6. In the current case, this corresponds to 
/i ~ 2.5 GeV. 
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FIG. 10: Ft versus (pa)^ with different masses. It has no significant mass dependence 
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The renormalized operators are defined as 



ZoO 



(52) 



The fact that the bare operator is independent of the renormalization scale //^ gives the 
renormalization group (RG) equation, 



2 d Q 



where 



7o 



2 '-^rcn • 



2/i2 dZo 



Zo diJ? 
is the anomalous dimension. 

The solution of Zoifi"^) can be written in the following form 

Co ifi'' 



(53) 
(54) 

(55) 



Zo (/i') 



^ > 7 fi/' 



Co if^") 

Expanding the anomalous dimension in the coupling constant 



(56) 



7o 



(57) 



and considering the running of in QCD f3 function (3{as) perturbatively 



An 



i=0 



(58) 



where are the coefficients of the QCD f3 function (3{as), one can solve for the coefficient 
functions Co perturbatively. 

The four loop solution j48| of the coefficient function in Eq. (jS^ is (we have suppressed 
the subscripts for the specific operator O) 



1+1 —. I (7i - Pi7o) 



TT 



An 



+ 



1 f as{ij) 
An 

An 



ill - f3iloT + 72 + /5i%o - /5i7i - /527o 



' 1 1 2 

^(fi - Pilof + 7^(71 - A7o)(72 + A To - Pill - /^27o) 
2 



+ ^ ( 73 - Pi^lo + 2/?i/327o - Psio + Pi^li - P2I1 - A72) } 



(59) 
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where 



_ ^ 7^ ^ A 



(60) 



The following tables give the anomalous dimensions_ 7_^^^ for Z^, Zm and Zt in the 



2G 



RI/MOM scheme for the quenched approximation 

Zs = Zp = 1/Zm, SO that C5(/i2) = Cpifx"^) = Note that, in Refs. 

the definition of Z's is the inverse of our definition in Eq. 



49|. In the case of chiral fermions, 



TABLE III: Quenched Anomalous Dimensions 7^ 



(0 



.(0) 



7 



(1) 



7 



(2) 



.(3) 



44.6667 



2177.0737 



130760.2969 



TABLE IV: Quenched Zm Anomalous Dimensions 'jm 



j{0) ^(1) ^(2) ^(3) 

8 252 11769.5469 557837.9375 



TABLE V: Quenched Zt Anomalous Dimensions 7-^ 



^(0) 


^(1) 


^(2) 


2.66667 


80.4444 


3268.2996 



The coupling constant itself is running with respect to /i. The four loop formula is given 
by 

^ 1 A In In (/xVA^c^) 

^P^H^J^ ^""^ (/^VAIcd) - P! In In (/iVA^cD) 

\n' In (/iVA^cD) " (2/3? - S^iP,) In In (/^VA^ci.) 



+ 2 
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The QCD /3-function is scheme independent only up to two loops. The additional terms 
of the expansion have been computed in the MS scheme in Ref. j5l| . 

In this work, the value of as was calculated at four loops using a lattice value of Aqcd 
taken from Ref. [^^1 as 

Aqcd = 238 ± 19 MeV . (62) 
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FIG. 11: Ta (labelled as "bare") and the scale invariant T^^ versus (pa)^ in the chiral limit. They 
coincide near (pa)^ = 4.1. The later is almost (pa)^ independent after (pa)^ > 2.4, the slope versus 
(pa)2 is about 0.001. 



Both Za and Zy are scale independent, but this is not the case for Z^. The scale invariant 
(SI) vertex for the axial and vector current is defined by removing the renormalization group 
running of Z^ as 

iM") = ^A,v iiapY) /C^ {{apf) , (63) 

where is defined in Eq. ()59j) with the anomalous dimension coefficients from Table IIIII 
We normalize C^{{^af' = 4.15) = 1, which corresponds to fi = 2.0 GeV, in order to compare 
with the {pay dependence of y ((ap)^). 

Fig, im shows both Ta {{apY, ma = 0) and F^^ ((op)^, rna = 0) as a function of (pa)^. By 
comparing them, we see that the renormalization group running due to Z^ is not appreciable 
for (pa)"^ > 3, but it does tend to make the SI data flatter as a function of {pay. The 
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FIG. 12: The same as Fig. JIT} for Tf^ and r™'''^^ The slope of SI versus {paf is about 0.008 
beyond (pa)"^ > 4.0. 

remaining scale dependence of {{ap)"^, ma = 0) is very small. A plausible explanation for 



it is an (ap) error [9|. Fitting the remaining scale dependence to the form 

T'nM')=^A+c{ap)\ 



(64) 



for a range of momenta that is chosen to be "above" the region where the condensate effects 
are important, one can obtain the scale invariant F^ which is denoted as F^^ 

When a linear fit of the SI data versus (ap)^ is performed for 2.4 < (ap)"^ < 5.7, the 
gradient is ~ 0.001. In the ideal case, the gradient should be zero. This small value is thus 
interpreted as an O(a^) error. This shows that the (pa)"^ error of the ratio of and is 
small, but we don't know their individual (pa)^ error separately. It appears that the (pa)^ 
error in as defined from the quark propagator in Eq. p9|) is as large as ~ 10% at p = 2 
GeV 0, Ql- However, this relatively large (pa)^ error in Z^^ must be cancelled to a large 
extent by that of Za, resulting in a small (pa)"^ error in F^^ (only 0.4% at {pa^ = 4.1). As 
will be discussed in the next sub-section, we will use Za determined from the chiral Ward 
identity to obtain Z^. Since Za, in this case, is determined from the pion state at rest, it is 
at a momentum scale of Aqcd- Thus, it should have small (pa)^ error. Using this value of 
Za and F^^ {{ap)^,ma = 0) to obtain and thus Z^\2 GeV) should give a (pa)^ error of 
~ 0.4% which we shall consider as the systematic error in 0{a?). 
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FIG. 13: The same as Fig. (jllj) for Tt and Tj^. The later is almost (pa)^ independent after 
(pa)^ > 2.0, the slope versus (pa)^ is about -0.0004. 



In the case of T^^^ = Z^/Zs, both Z^ and Zs run with /x^. Fig. shows F™^ and the 

sul 
S 



corresponding scale invariant (SI) vertex, F^" ' , after three loop running. We see that F""^'^' 



is much flatter than F^"*^ for {pa)"^ > 3.0. The linear fit of the SI data versus (pa)^ in the 
range of 4.0 < (ap)^ < 5.7, gives a gradient of 0.008(3). This is an order of magnitude larger 
than that of the axial vector (and that of the tensor current below). This relatively larger 
gradient could be due to the systematic uncertainty in subtracting the chiral divergence in 
F5 or the mismatch of the (pa)^ errors between Z^^ and or both. 

The SI result of Ft which comes from the three loop running of Zt j3| and four loop 
running of Z^ is plotted in Fig. along with Ft- The linear fit to the SI data in the range 
of 2.0 < {apf < 5.7 gives a gradient of -0.0004(2). 

It is worthwhile pointing out that comparing to the Domain Wall fermion case on a 



16^ X 32 X 16 lattice with Wilson gauge action at /3 = 6 



.0 3, 



we find that the remaining 



scale dependence as measured by the gradient in (pa)^ is comparable for the F|^ case. But 
the gradient for F^^F^^ is one/two orders of magnitude smaller than that in the Domain 
Wall fermion case. Since the Domain Wall fermion with finite 5th dimension and the rational 
polynomial approximation of the sign function that we adopt in the present work are two 
different approximations of the same overlap fermion in 4 dimensions j^l, the (pa)^ errors 
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in the scale invariant vertices serve as a measure of the O(a^) errors of the approximation. 



A. Determining from Ward identity 

Before applying the renormalization group running (an inverse operation of Eq. (j63j) ) to 
match to the MS scheme at certain scale, we need to input to Eq. (f^Hj) in order to de- 
termine other renormalization constants from the respective vertex functions. As explained 
in Sec. Ill Bl we prefer using Zp^ from the chiral Ward identity to determine Z^ than directly 
obtaining it from the quark propagator. This is partly due to the fact that there is ambiguity 
in the lattice definition of momentum lof in Eq. (11911 . Furthermore, it is shown in the study 

n n 

with Domain Wall fermion [£| and an earlier study of the overlap fermion that the (pa)^ 
errors in the scale invariant Z"^ are quite large. 

The renormalization constant Zj^ for the axial current = — D/2p)^)ip can 
be obtained directly through the axial Ward identity 

ZaO^A^ = 2Z^mZpP, (65) 

where P = tpiv y^jl — D /2p)^)ip is the pseudo-scalar density. For the case of the overlap 



fermion Q,Q9|, |2l|, Z^ = Zg^ and Zs = Zp. Thus, Z^ and Zp cancel in Eq. and one 
can determine Z^ to O(a^) non-perturbatively from the axial Ward identity using the bare 
mass m and bare operator P. To obtain Za, we shall consider the on-shell matrix elements 
between the vacuum and the zero-momentum pion state for the axial Ward identity 

ZA{0\d^A^\7T{p = 0)) = 2m(0|P|7r(p = 0)), (66) 

where the matrix elements can be obtained from the zero-momentum correlators 

Gd,A,piP = 0,t) = {J2diAiix)P{Q)) 

Gppip = 0,t) = (X]P(x)P(0)). (67) 

X 

The non-perturbative Za is then 

m-.0,t^oo Ga^A4P{P = 0,t) 

Given that the time derivative itself in Gd4^A4p{'P = 0,t) invokes an O(a^) error, it would 
be better to adopt a definition for Za which is devoid of this superfluous O(a^) error. This 
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FIG. 14: ZA{ma) vs quark mass ma. 

can be achieved by noticing that, at large t where the pion state dominates the propagator 
GdiAiPip = 0,t), one can effectively make the substitution 

Go,aAp = 0,t)^-^^m^GAAP = 0,t). (69) 



Consequently, Eq. (jHI^ becomes 

m^O,i^oo rriT^GAiPip = 0, t) 

Since the Gpp and G'a4P correlators are calculated at finite ma, we shall define the 
renormalization factor 

7 f ^ r '2mGpp{p = 0,t) 

ZA(ma) = lim — -— -, (71) 

where the massless limit would give Za- We plot the results of Z^(ma) from Eq. ()71|) in 
Fig.^S In view of the fact that there is no 0(a) error with the overlap fermion, Zyi(ma) could 
have terms like mAgc^a^ and m?a?, but could also have terms like m/ Kqcd and (m/AornY- 
Thus, we parameterize it with the form which is linear and quadratic in m 

Zyi(ma) = Za{1 + bAm + CArn^)- (72) 

As noticed before j2l|, it is conspicuously fiat as a function of ma suggesting that the O(a^) 
error from the action and the axial-vector operator is small. From the fitting to the form in 
Eq. (1721) for the range of ma from 0.014 to 0.6, we find that Za = 1.849(4), Ba = -0.347(16) 
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(in units of Kqcdo^ with Aqcd = 0.238 GeV), and ca = 0.317(21) (in units of a^) with 
xVdof = 0.54. 

We observe that Za is determined to the precision of 0.2% in statistical error. It is thus 



more desirable j9|| to use Za from the Ward identity and Ta to determine all the other 
renormalization constants. 

From the renormalization condition Eq. ()22|). we finally obtain Z^^ 



= Za^,. (74) 



Zf = ZaT^I (73) 
and the other SI renormalization constants from Eq, 

pSI 

o 

Since we normalize the scale invariant vertex functions to the RI scheme at (pa)^ = 4.1, 
the renormalization constant determined in Eq. ()74|) is just Zq^{^ = 2 GeV). 

B. Matching to MS scheme 

In order to confront experiments, one frequently likes to quote the final results in the 
MS scheme at certain scale. For light hadrons, the popular scale is 2 GeV. To obtain 
the renormalization constants in the MS scheme at 2 GeV, one can use the perturbatively 
computed coefficient functions in Eq. in the MS scheme to evolve the scale invariant 
renormalization constant to the targeted scale, i.e. 

Z^if,) = CW{1^)Z'^. (75) 

Alternatively, one can avoid the step of going to the scale invariant quantity and, instead, 
match directly from the Rl-scheme to the MS scheme at the same scale. The perturbative 
expansion of the ratio Z^^/Z^^ to two-loop order are given by the finite coefficients in the 
perturbative expansion of Z^^ 

= 1^ = 1 + ^{Z^^)^' + i^nz^^r + •••• (76) 
The numerical values of the matching coefficients, Z^^^ , Z^ and Z^^^ in Eq. (j76|) used for 



Z^ , Zm and Zt have been calculated in Ref. j48[ and Ref. [5J|, and are collected in Table EH 
The results for Za, Zy, Zp, Zs, Zt and Z^ in the RI and MS scheme at 2 GeV (for 
Zpi Zg^ Zt and Z^) are listed in Table IVlll 
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TABLE VI: Quenched RI to MS matching coefficients 



Z(o) Z};' Z^q"> Z^q'> R at 2 GeV 

0.0000 -25.4642 -1489.9805 0.98706 

Zm -5.3333 -149.0402 -5598.9526 0.85127 

Zt 0.0000 -46.6654 -2067.9753 0.97909 



TABLE VIL RenormaUzation constants Z in RI, and MS schemes at /i = 2 GeV. The renormal- 
ization constant in the SI scheme is normaUzed to be the same as that in the RI scheme at 2 GeV. 
These renormahzation constants are obtained from the lattice with a = 0.200 fm. 



z 


RI at 2 GeV 


MS at 2 GeV 


Za 


1.853(9) 


1.853(9) 


Zy 


1.846(8) 


1.846(8) 


Zp 


1.571(15) 


1.845(17) 


Zs 


1.567(13) 


1.841(15) 


Z'x 


1.966(6) 


1.925(6) 


z^ 


2.307(18) 


2.277(18) 



VI. FINITE m DEPENDENCE OF THE RENORMALIZATION FACTORS 



As discussed in Sec. IIIIl a non-perturbative renormalization of the heavy-light axial cur- 
rent via the cliiral Ward identity and the unequal mass Gell-Mann-Oakes-Renner relation is 
possible with the overlap fermion at finite ma This offers an opportunity to calculate 
heavy-light decay constants and transition matrix elements without being subjected to the 
uncertainty of the perturbative calculation of the renormalization constants. In this section, 
we shall examine the finite m behavior of the renormalization factors. This is useful for the 
future study of heavy-light decays and transitions. 

Similar to the renormalization factor Zyi(ma) defined in Eqs. (1711) and (f7^ for the axial 
current, we shall parameterize the renormalization factor for the other operators by 

Zoirna) = ZoifJ^a, g{ci)){l + bom + Cpm'^), (77) 

where the linear m term includes terms like mKqcDO^'^ and m/ Kqcd and the quadratic m 

32 



term includes terms like m^a^ and {m/ KqcdY ■ 

To obtain the renormalization factors for O = V,S,P,T and ip, we shall adopt the same 
renormalization condition Eq. (jl4|) for finite m. We shall use the quark propagator to obtain 
the finite m dependence for Z^{ma) from Eq. (jl9|) . except with = Zr^{ma = 0) normalized 
from Za via the renormalization condition, Eq. (j22j) . We then follow the above procedure 
in the preceding sections to obtain ZAifna), Zv{ma), Zsima), Zp{ma), and ZTima). The 
Z^^{nia, /i = 2 GeV) for Z^{ma) in the MS scheme at /i = 2 GeV is plotted in Fig. 1201 The 
corresponding results on ^^(ma), Zy^ma), Z^^{ma,fi = 2 GeV), Z^^ima^n = 2 GeV) and 



ma. 



/i = 2 GeV) are plotted in Fig. Cg CEl CHI El and d respectively. 



We observe that these curves are all rather flat, with less than 3% deviation for ma as 
large as 0.6. We make a correlated fit of the finite m behavior with the form in Eq. ()77j) . 
The coefficients bo and cq are listed in Table IVIIII The fitted curves are drawn in Figs. El 
El El El and El as solid lines. For comparison, we also plot the uncorrelated fits as dashed 
lines which have much smaller x^/dof than those of the correlated fits. The flatness in 
ma suggests that the OiimaY) dependence is weak. This is consistent with the findings 
for the vr and p masses l^f and their dispersion relations 12|. The observed 0{ma^) and 



0{m a ) dependence are much weaker than those found with the chirally improved Dirac 
operator [10|. We don't know the finite m behavior of the domain- wall fermion from the 
study of the quark bilinear operators ^. The study of Z^(ma) and Zy{ma) from the nucleon 
matrix elements found much stronger 0{ma'^) and 0{m'^a'^) dependence than observed here. 

We should point out that in Sec IV A[ ZAijna) was defined in Eqs. (ffT|l and (f7^ with the 
implicit assumption that Zp{ma) is not very different from Zs{ma) for the available range 
of ma. Even though this does not in principle affect the extraction of the renormalization 
constant Za at the massless limit of Z^(ma), where Zp is equal to Zs, the extrapolation 
could have a large uncertainty if Zp{ma) is very different from Zs{ma) for the range of ma 
which is not close to the chiral limit. Thus, it is gratifying to see in Figs. ITTI and ITSl that 
the ma behaviors of Zp{ma) and Zs{ma) (in the MS scheme at 2 GeV in this case) are 
very similar. They differ less than 3% for ma < 0.6. The upshot is that we should take the 
ma behavior in Zyi(ma) in Fig. El as the correct one, instead of that in Fig. ^^as defined 
in Eq. ()71|) . Any attempt to correct for the negligence of the ma dependence in Zp{ma) 
and Zs{ma) in Eq. (f7T|) and compare with that in ZAifna) in Fig. El is expected to reflect 
different finite ma corrections under different renormalization conditions. 
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TABLE VIII: Mass dependence of the renormalization factors as defined in Eq. ()77() . bo is in units 
of Aqcdo^ with Aqcd = 0.238 GeV. cq is in units of a^. 



Z^S (2 GeV) 



bo 



CO 



xVdof 



A 
V 
P 
S 
T 



1.853(9) 
1.846(8) 
1.845(20) 
1.841(15) 
1.925(6) 
2.277(18) 



0.029(16) 
-0.0055(33) 
-0.015(12) 
0.22(19) 
-0.012(8) 
-0.001(2) 



-0.003(2) 
-0.013(7) 
0.017(4) 
0.004(2) 
-0.015(3) 
0.002(1) 



1.3 
2.2 
0.3 
0.7 
1.6 
0.8 




FIG. 15: Renormahzation factor ZA_{ma) against quark mass ma. The soHd line is the correlated 
fit and the dashed hne is the uncorrelated fit. 

VII. SUMMARY AND OUTLOOK 



In this work, we performed a non-perturbative renormalization calculation of the com- 
posite quark bilinear operators with the overlap fermion in the regularization-independent 
scheme from the quark vertex function with high virtuality. The renormalization group 
running of the renormalization constants were calculated to obtain the scale invariant (SI) 
renormalization constants and also matched to the MS scheme at = 2 GeV. The scale 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 

ma 



FIG. 16: The same as in Fig. ^Jfor Zy(ma) 



2.05 




0.6 0.7 



- MS, 



FIG. 17: The same as in Fig. El for the renormahzation factor Zp {ma) in the MS scheme at 
fi = 2 GeV. 



invariant Za from the Ward identity and the axial vertex function were used to eliminate 
the wavefunction renormahzation Z^ and determine the renormahzation constants from 
other vertex functions. Since Z^, as obtained from the Ward identity, has a very small 
error (~ 0.2%), it is more desirable to use it to determine Z^ instead of using the quark 
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2.05 
2 

1.95 



T r 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



FIG. 18: The same as in Fig. Elfoi' Zs^^ima 




FIG. 19: The same as in Fig.HZIfor ZT^^{ma). 

propagator to determine it. The latter can introduce an error as large as ~ 10 - 20% j^. 

After subtracting the quenched chiral log divergences in the vertex functions Tp and Ts 
due to the presence of the pseudoscalar meson, the expected relation Zs = Zp due to chiral 
symmetry holds to high precision (~ 1%) for a large range of (pa)^ with (pa)^ > 3. The same 
is true for the relation Za = Zy. The resultant check on the chiral symmetry relations are 
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2.5 
2.45 
2.4 



2.35 

<v 
O 

<^ 2.3 

II 

2,25 



2.2 
2.15 
2.1 



T 1 r 



J I L 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 



FIG. 20: The same as in Fig.HHfor Z^^^{ma). 

comparable to those of the domain-wall fermion P] and somewhat better than the chirally 
improved fermion where it is found that Z^/Zy ~ 1-03 and Zp/Zs ~ 0.95 for their 
smallest lattice spacing at a = 0.078 fm. 

We studied the finite m behavior in the renormalization factors of these composite oper- 
ators. This is useful if one wants to use the same overlap Dirac operator for both the light 
and heavy quarks. With present-day computers, it is not practical to reach a lattice spacing 
such that ma <^ 1 for the charm and bottom quarks. As such, one would like to have a 
Dirac operator which has chiral symmetry and, at the same time, has small 0{ma?) and 
0{m?'a^) errors. It is suggested |l2| that the overlap fermion with its effective quark propa- 
gator having the continuum form might be suitable for this purpose. Indeed the 0{ma^) and 
0{m?a^) errors in the dispersion relation are shown to be small [17J. In this study, we find 
that the finite m dependence is quite gentle in Z^{ma), Zyi{ma), Zv{ma), Zp{ma), Zs{ma), 
and ZxiTna). For ma as large as 0.6, the deviations are generally less than 3.5% in these 
renormalization factors we studied. 

Since the lattice community is geared to carrying out large scale djTiamical fermion 
calculations with chiral fermions, it is worth studying the (pa)^ errors in the scale invariant 
vertices to gauge the O(a^) errors and also the mass dependence of the renormalization 
factors with the overlap fermion. 
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